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We present a modelling approach for diffusion in a complex medium characterized by a random 
length scale. The resulting stochastic process shows subdiffusion with a behavior in qualitative 
agreement with single particle tracking experiments in living cells, such as ergodicity breaking, p- 
variation and aging. In particular, this approach recapitulates characteristic features previously 
described in part by the fractional Brownian motion and in part by the continuous-time random 
walk. Moreover, for a proper distribution of the length scale, a single parameter controls the ergodic- 
to-nonergodic transition and, remarkably, also drives the transition of the diffusion equation of the 
process from non-fractional to fractional, thus demonstrating that fractional kinetics emerges from 
ergodicity breaking. 
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I. INTRODUCTION 

Many processes in life sciences, soft condensed matter, 
geology, and ecology show a diffusive behaviour that can¬ 
not be modelled by classical methods. These phenomena 
are generally labeled with the term anomalous diffusion 
in order to distinguish them from the normal diffusion, 
where the adjective normal has the double aim of high¬ 
lighting that: i) a Gaussian-based process is considered; 
and ii) that it is a usual diffusion process with a linear 
growth in time of the particle displacement variance. The 
observation in nature of anomalous diffusion has been 
definitively established experimentally and several theo¬ 
retical models have been proposed for the interpretation 
of such phenomenon [1-3]. Among these theoretical ef¬ 
forts, the fractional calculus has emerged to be a suc¬ 
cessful tool for modelling a class of anomalous diffusion 
processes [4, 5]. For this reason, anomalous diffusion gov¬ 
erned by equations built on fractional derivatives is often 
also referred to as fractional diffusion. Several stochas¬ 
tic approaches have been proposed in the literature to 
reproduce fractional kinetics [6-10]. 

In the last decades, advances in fluorescence-based 
techniques such as single-particle-tracking (SPT) have 
allowed to precise characterization of the diffusion of 
molecules in biological systems [11]. In particular, the 
recording of long single-molecule trajectories has revealed 
that the occurrence of anomalous diffusion of some cellu¬ 
lar components in living cells is associated with ergodic¬ 
ity breaking (EB) [12-16], i.e. the nonequivalence of time 
and ensemble averages [2, 3]. Often EB and anomalous 
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diffusion are concomitant with aging, i.e. the dependence 
of statistical quantities on the measurement time [17]. 

Besides the fundamental interest of nonergodic pro¬ 
cesses in statistical mechanics and its still unclear impli¬ 
cations in cell biology, the occurrence of EB further em¬ 
bodies a valuable criterium for the selection of the under¬ 
lying diffusive stochastic process. In this respect, com¬ 
parative studies - involving the fractional Brownian mo¬ 
tion (fBm) [18, 19], the fractional Langevin equation [18] 
and the continues time random walk (GTRW) [19-21] - 
have been conducted in order to determine which type of 
motion could possibly cause nonergodic anomalous diffu¬ 
sion. Among the mentioned theoretical frameworks, the 
fBm and the fractional Langevin motion are ergodic, with 
fBm displaying EB only in the ballistic limit [18]. On the 
other hand, the GTRW is nonergodic [20, 21] with the EB 
stemming from the nonstationary nature of the process 
when the distribution of waiting times has a power-law 
tail [20]. Eor this reason, the GTRW has been extensively 
used to model the occurrence of nonergodic diffusion and 
the waiting times have been associated to immobilization 
events caused by biochemical interactions [2, 3]. 

However, due to the lack of nonergodic models alterna¬ 
tive to the GTRW, the use of EB as a criterion to select 
the dynamic process has shown some limitations. An 
example is provided by the seminal work of Golding and 
Gox [12]. In this case, although the presence of EB favors 
GTRW as the model underlying the dynamics of RNA in 
cellular cytoplasm, a moments-based criterion called p- 
variation [19] seems to indicate a diffusion compatible 
with fBm. Similarly, other experiments also showed the 
simultaneous occurrence of EB and nonlinear scaling of 
the time-averaged mean-squared displacement, making 
necessary to hypothesise the coexistence of GTRW with 
other processes in order to theoretically model the ob¬ 
served features [13-15]. 
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In this paper, we provide a general framework in which 
EB emerges as a consequence of the heterogeneity (or 
randomness) of the system. The heterogeneity is de¬ 
scribed by the random nature of a characteristic property 
of the medium, such as a length scale depending on a 
single parameter (3. Simple examples of this behavior are 
provided by a population of particles, each of them diffus¬ 
ing in a Brownian fashion but with a broad distribution 
of diffusion coefficients However, our conclusions do 
not depend on the type of motion performed by the par¬ 
ticles. We also show that for any nontrivial choice of the 
distribution ip, the parameter /3 continuously drives the 
transition from ergodic to a nonergodic process. Notably, 
a fractional kinetics straightforwardly emerges from EB, 
and thus allows us to associate nonergodicity to a frac¬ 
tional equation. 

For its generality, our approach constitutes a flexible 
tool to interpret the occurrence of EB in random media 
and in living cells without involving CTRW and subor¬ 
dination. From the biophysical point of view, it implies 
that EB can be generated by heterogeneity in the diffu¬ 
sion, without the need of particle trapping. In partic¬ 
ular, we discuss how our model can resolve the contro¬ 
versy on the interpretation of Golding and Cox experi¬ 
ments [12, 19, 21, 22] by considering the fBm in a het¬ 
erogeneous medium. Such a model allows one to simulta¬ 
neously obtain the apparently contradictory features ob¬ 
served in Golding and Cox experiments, i.e. the mono- 
tonically increasing p-variation test typical of the fBm 
together with the EB parameter of the CTRW. 

Finally, we show that our formulation can be further 
generalized by considering a nonstationary length scale 
if} = and thus including the occurrence of aging. 


II. ERGODICITY BREAKING FROM 
DIFFUSION IN A RANDOM MEDIUM 

In our model, we consider a stochastic process defined 
as: 


X{t) = ip Xgen{t) , 

describing a population of particles diffusing according 
to a generic ergodic Gaussian process Xgenii) in a com¬ 
plex random medium. The medium properties are inde¬ 
pendent of the diffusing particles and its randomness is 
described by a random characteristic quantity - such as 
a length scale ip - with distribution depending on the 
parameter (3. The role of (3 thus consists in tuning the 
degree of randomness of the medium by modulating the 
distribution of the length scale. For for the case in which 
Xgen{t) represents a random walk, ip corresponds to a 
distribution of diffusion coefficients. 

Although the following conclusions hold for every er¬ 
godic Gaussian stochastic process, for the sake of simplic¬ 
ity from now on we will consider the fBm X^it), an er¬ 
godic non-Markovian Gaussian process characterized by 


the covariance matrix: 


^Hit,s)=t^^ + , ( 1 ) 


where 0 < i? < 1 is the Hurst exponent, and the variance 
results to be (Xj^) = . 

Therefore, we investigate the following diffusion pro¬ 
cess A(t) in a random medium: 

X{t)=ipXH{t). (2) 

In order to study the dynamics of the process, we 
first consider the time-averaged mean-square displace¬ 
ment [18, 20, 21] 

where A is the timelag and T the measurement time. 
The time-averaged mean-square displacement describes 
the time dependence of the second moment of the parti¬ 
cle’s position and it is often used to classify the diflh^ion 
mode. For the pure Brownian motion (2iJ = 1), 6^iT) 
shows a linear growth with A, whereas for the fBm shows 
a power-law behaviour ~ A^^, i.e. anomalous diffusion. 
The effect of the random length scale is preserved in the 
calculation of (5^(T). For the particular case 2H = 1 in 
which the process Xnit) in (2) corresponds to the pure 
Brownian motion, the random length scale is propor¬ 
tional to the diffusion coefficient. Gonsequently, as shown 
in Fig. 1, time averages such as 5^(T) remain random 
variables and thus irreproducible [23], causing ergodicity 
breaking (EB). This effect can be estimated through the 
calculation of the EB parameter Eb{T) [18, 21]. Let (•) 
represent the ensemble averaging, then 


Ab(T) 


mT)?) ^ 


( 4 ) 


is calculated in the large T limit and tends to 0 when 
the process is ergodic [18]. 


With a fixed and non-random length scale, e.g. ip = 1, 
for the stochastic process A(t) defined in (2), we obtain 
[18] 




\t) 


T—>-oo 


t 0. 


( 5 ) 


In contrast, if - 
that 


is a random variable, for X{t) it holds 




i^t) 


& 








( 6 ) 

The condition (£^) > (£^)^ is met in general for any 
distribution as a consequence of the inequality K > -1-1 

[24], where K and S are the kurtosis and the skewness 
respectively, and in particular for any unilateral non¬ 
increasing density it holds K >9/5 from the Gauss- 
Winckler inequality [24]. The limiting case (£|) = (^|)^ 
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is met when the distribution of the length scale is the 
Bernoulli distribution with equal success probability for 
values 0 and 1 or it is the Dirac-delta distribution 
6 {if} — 1); therefore the process is nonergodic for every 
nontrivial choice of Although these conclusions might 
look somehow trivial, they show how a complex medium - 
through a random distribution of the length scale - might 
produce nonergodic behavior into an ergodic Gaussian 
stochastic process, including the pure Brownian motion, 
only by introducing heterogeneity [20]. 


III. ERGODICITY BREAKING AND THE 
FRACTIONAL KINETICS 

In the previous section we have shown that, since 
is an independent random variable, EB can occur as the 
sole consequence of the randomness of the medium in 
which diffusion takes place (6) and independently of the 
chosen ergodic Gaussian stochastic process. 

In the following, we will focus our attention on the 
stochastic process X{t) as defined in (2). This process 


has already been studied in a specific characterization 
named generalized grey Brownian motion (ggBm) [8, 25, 
26]. As a matter of fact, the ggBm trajectory is 

obtained by setting ip = i.e. 

= \/^Xnit), (7) 

where the positive random variable Ap is distributed 
according to the one-side M-Wright/Mainardi function 
Mp{X), with A > 0 and 0 < /3 < 1, defined as [27, 28] 

The case of a non-random length scale, i.e. Ap = I, is 
straightforwardly recovered in the limit /3 —>■ I since it 
holds Ml (A) = 6{X — I). The ggBm is a rather general 
model and includes as special cases the Brownian motion 
(/3 = 2H = 1), the fBm {/3 = I) and the grey Brownian 
motion (/? = 2H). 

It is well known that the probability density function 
of Xp^nit) is [8] 




1 

Ay(27rA)”det jh 





Mp{X) dX, 


(9) 


where x = (xi,...,a;„) and jh = i,j = 

1,..., n, is the covariance matrix of the fBm defined in 
(1). Therefore, by the Mellin transform of Mp{X) [29], 
i.e. /“ X^-^Mp{X) dX = r[l + (s - l)]/r[I + - 1)], 

with s > 0, the covariance matrix of the ggBm can be 
obtained as [8, 26] 

1/3,nit, s) = -\t- sp") • (10) 

The one-point one-time density function can be derived 
from (9) and becomes 

= ( 12 ) 

where it emerges that the shape of probability density 
function of displacements is affected by the medium, here 
represented by Mp{X). In terms of the H-function the 
density function V{x]t) reads [30, 31] 


10 






- ; (I-/3/2,/3/2) 

(0,1) ; - J ’ 

(13) 


and the asymptotic decay is Mp/ 2 {\x\ —>■ oo) ~ 
|a;|f(/3-i)e-f'|iiib, with b = and c = 2 ^ [31, 


32]. From (10) the variance turns out to be 


(^Ih) - 


J.2H 


(14) 


showing that the presence of the medium does not af¬ 
fect the power law growth of the particle displacement 
variance over time. It is noteworthy to observe that the 
ggBm shows both subdiffusion, 0 < M < 1/2, and su¬ 
perdiffusion, 1/2 < < 1. Moreover, a remarkable case 

is represented by M = 1/2 in which the particle displace¬ 
ment variance results to be linear in time, see (14), but 
the density function is not Gaussian according to (12). 
The Gaussian density is obtained from (12) as a special 
case when /3 = 1. 

The evolution equation for V{x;t) is given by 


dt (3 21 //^ Q^2 ’ 


(15) 


where is the Erdelyi-Kober fractional derivative 

with respect to t and then the process is also referred to 
as Erdelyi-Kober fractional diffusion [33]. Special cases 
of Eq. (15) are: the classical diffusion (/? = 2H = 1), 
the fBm master equation (/3 = 1) and the time-fractional 
diffusion equation {j3 = 2H). A similar approach can 
be developed in the framework of the space-time frac¬ 
tional diffusion equation, which includes all its special 
cases [34]. 
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We would like to remark that the fractional kinetics, 
i.e. /3 ^ 1, emerges directly from the EB due to the 
randomness of since M^^i(A) 7 ^ 5{\ — 1). 

Moreover, the fractional order related to (3 can be ex¬ 
perimentally computed by means of the long-time limit 
of the EB parameter. In fact, for large T, from ( 6 ) and 
(.jj = the EB parameter then becomes 


E 


(ggBm) 

B 






r(2/3) 


(16) 


where again the Mellin transform of M^^X) [29] has been 
used to compute (A^) and (A^). 

In summary, the existence of a random length scale 
turns an ergodic process into a nonergodic one without 
the need to introduce an alternative stochastic process. 
When this transition occurs continuously with respect 
to a parameter /3, the distribution of the length scale 
can be related to the M-Wright/Mainardi function and 
the resulting stochastic process is driven by a fractional 
diffusion equation. 

Therefore, the present formulation provides a founda¬ 
tion of fractional kinetics on the basis of the appearance 
of the EB. In other words, fractional kinetics can be con¬ 
sidered as stemming from the EB due to the heterogene¬ 
ity of the medium in which the diffusion takes place. In 
order to support this physical foundation argument, we 
remark that from the proposed ggBm (7) the evolution of 
the particle density function is governed by a fractional 
diffusion equation also in the special case H = 1/2, see 
(15), with Xnit) performing the classical Brownian mo¬ 
tion and the particles displaying a variance with a linear 
growth in time (14). 


IV. RELATION WITH EXPERIMENTS 

Advances in biophysical techniques, such as SPT, 
have allowed researchers to detail the motion of sin¬ 
gle molecules and have revealed very complex diffusion 
patterns in living-cells [ 11 ]. In particular, the analysis 
of these experiments has shown that several biological 
systems display nonergodic behavior as a consequence 
of interactions occurring in heterogeneous cellular en¬ 
vironments [12-16]. Such nonergodic behavior has of¬ 
ten been connected with the occurrence of anomalous 
(sub)diffusion. The occurrence of EB has been mainly 
identified through the nonequivalence of time and en¬ 
semble averages and by the calculation of the EB param¬ 
eter (4) [18, 21]. Owing to the importance of molecular 
transport for the cellular function, theoretical efforts have 
been devoted to understand the physical mechanism be¬ 
hind EB in biology. Several stochastic models presenting 
nonstationary (and thus nonergodic) (sub) diffusion have 
been proposed [3]. Among these models, the most popu¬ 
lar has definitively been the CTRW [19-21, 35, 36] which 
has been extensively used to model nonergodic subdiffu¬ 
sion in living cells [13-15]. The CTRW has allowed asso¬ 
ciation of the nonergodic behavior with the occurrence of 


particle immobilization with a heavy-tailed distribution 
of trapping times [37]. 

However, among the experimental evidences of EB in 
biological systems, not all the observed features could be 
directly addressed within the framework of CTRW alone. 
For example. Refs. [12, 14, 15] showed subdiffusive scaling 
of the time-averaged mean-square displacement obtained 
for single trajectories, making necessary the postulation 
of the coexistence of CTRW with other sources of subd¬ 
iffusion, i.e. the fBm [15] or a fractal processes [14], in 
order to properly interpret the results. In addition, some 
experiments did not show the occurrence of inherent fea¬ 
tures of CTRW, such as aging [38] or immobilization [16]. 

In order to determine the physical scenario behind the 
subdiffusive EB, a number of diagnostic tools have been 
proposed [39]. Among these, a valid criterion for selection 
of stochastic processes is represented by the so-called p- 
variation test [19]. The test is based on the calculation 
of the quantity 


V^P\t) = \im , (17) 


where for t € [0, T] 


2^-1 


j=o 


X 


U + 1 )^ 

2 " 




(18) 

with a Ab = minja, b}, and allows the CTRW-like mod¬ 
els and the fBm to be distinguished, even on the single 
trajectory level [39]. 

In spite of the efforts in developing tests and meth¬ 
ods to distinguish between different stochastic models, 
contradictory indications still prevent the unambiguous 
determination of the physical mechanism behind EB in 
biological samples. An example is provided by what is 
probably the first evidence of EB in living cells, i.e. the 
experiments describing the motion of individual mRNA 
molecules inside living E. coli cells presented in the sem¬ 
inal paper by Golding and Cox [12]. In this case, in 
order to explain the occurrence of EB as evidenced by 
the large scattering of single-trajectory 6“^ curves and 
a non-zero EB parameter, the CTRW was proposed in 
Refs. [21, 22] to model this dataset. However, in or¬ 
der to account for the subdiffusive behavior of the time- 
averaged mean-square displacement, the authors of both 
works proposed the coexistence of CTRW with some 
degree of spatial confinement producing the power law 
behavior of 6^ [21, 22]. But the application of the p- 
variation test to the same dataset [ 12 ] showed that the 
subdiffusion is unlikely to originate from the CTRW, 
whereas the data are compatible with fBm [19]. 

In this scenario, the general stochastic process pre¬ 
sented in this work in ( 2 ) provides a plausible frame¬ 
work to describe the subdiffusive nonergodic behavior 
observed in Ref. [12]. The introduction of a random 
length scale associate to a random medium allows to de¬ 
scribe the complexity of the cytoskeletal environment and 
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reproduce the scatter of time-average mean-square dis¬ 
placement observed at the single trajectory level. This 
observation is quantitatively translated by the calcula¬ 
tion of the EB parameter. As a matter of fact, equa¬ 
tion (16) shows that the EB parameter of the specific 
process described in (7) is identical to the one obtained 
for a CTRW with a power law distribution of waiting 
times, i.e. ipir) oc and infinite average sojourn 

time [21, 22], independently of the ergodic Gaussian pro¬ 
cess used to model diffusion. In addition, the flexibility of 
our method allows us to choose the fBm to model single 
particle diffusion (7) and thus reproduce the subdiffusion 
in 6“^ and maintain the same p-variation behaviour of the 
^ggBm (^) = ^fBm (^) > while preserving the same 
degree of EB observed for CTRW-like models (Fig. 2). 

V. AGING 

An interesting feature emerging from some single¬ 
particle tracking experiments of cellular components [14- 
16] is the occurrence of aging, i.e. the dependence of sta¬ 
tistical quantities - such as the time and ensemble aver¬ 
aged mean-square displacement {d‘^{T)) - on the measure¬ 
ment time, as a consequence of the presence of nonsta- 
tionarity in the diffusive mechanism [17]. Besides living 
cells, aging has been observed for many complex systems 
such as blinking nanocrystals [40-42], spin glasses [43] 
and colloidal suspension [44]. Since aging can charac¬ 
terize long-term memory [45], it can be used as a sta¬ 
tistical indicator of complexity and thus exploited to dis¬ 
criminate among different modelling approaches [46, 47]. 
Furthermore, aging has been shown to be associated with 
weak ergodicity breaking [48, 49], i.e. a situation in which 
the time needed to explore a system phase space is infi¬ 
nite but the phase space can not be divided into mutually 
inaccessible regions [48]. 

Our theoretical formulation allows to reproduce aging 
by the extension to the case of a nonstationary random 
medium ip = £p{t) [50]. The stochastic process results 
to be defined as: 

Xa,0,H{t) = Xnit) , (19) 

where Ap and Xnit), with 0 < H < 1, have the same 
meaning as in equation (7). In this case, the increments 
of Xa^p^nit) are nonstationary, in contrast to the process 
defined in (7), which is recovered as a particular case for 
a = 0. The parameter a is constrained by the physical 
requirement that the process is diffusive, meaning that 
the particle displacement variance must grow in time. 
Since the variance of the process is given by 

= ( 20 ) 

the latter condition can be expressed as a > —2H. It 
can be shown [50] that the time- and ensemble-averaged 
mean-square displacement then is {6‘^{T)) ~ ^ 2 Hrpa 
(Fig. 3). 


It is interesting to note that our formulation shows 
properties that were not recapitulated by any of the mod¬ 
els for nonergodic diffusion previously presented in liter¬ 
ature [3]. First, the exponents controlling the power law 
behavior of d'^, ^ and {S'^{T)) depend on two pa¬ 

rameters, a and H. As such they can thus be indepen¬ 
dently tuned to reproduce any different scaling of the two 
curves, in contrast to the other models [3]. In particular, 
our model show that the time- and ensemble-averaged 
mean squared displacements can have marked different 
behavior, for example with one showing subdiffusivity 
while the other showing superdiffusivity. In addition, the 
aging can shows positive or negative exponent depend¬ 
ing on the relative magnitude of the exponents control¬ 
ling the growth of the time- and ensemble-averaged mean 
squared displacement (Fig. 4). 

Moreover, we highlight that the aging can be obtained 
even in the case in which the time-averaged mean squared 
displacement or the ensemble-averaged mean squared 
displacement (X^ p ^) show Brownian behaviour, i.e. 
when 2H = 1 or a = 1, respectively. It is interest¬ 
ing to note that in the case 2H = I we recover the 
same relationship between the exponent of the ensemble- 
averaged mean squared displacement (a -I- 1) and the 
time-ensemble-averaged mean squared displacement ob¬ 
tained for other models, such as the CTRW [20], the 
scaled Brownian motion [51], the quenched trap [52] and 
the patch model [53]. Moreover, the calculation of the 
EB parameter (4) for the process (19) shows that even in 
the presence of aging {a ^ 0) the value of the EB param¬ 
eter is identical to the one obtained for a CTRW with 
infinite average sojourn time and power law distribution 
of waiting times [50]. 


VI. CONCLUSIONS 

We have demostrated that an ergodic Gaussian process 
occurring in a heterogeneous medium characterized by a 
random length scale can be turned into nonergodic with¬ 
out altering the properties of the Gaussian process itself. 
We showed that for any nontrivial choice of the distribu¬ 
tion of the length scale, the transition from ergodicity to 
nonergodicity can be continuoulsy tuned by means of a 
parameter j3. In these cases, the distribution of the length 
scale can be related to the M-Wright/Mainardi function 
and the resulting stochastic process is controlled by a 
fractional diffusion equation. 

These conclusions are valid for any ergodic Gaussian 
process. Therefore, the generality of our formulation 
posits it as a flexible tool for the interpretation of het¬ 
erogeneous and/or nonergodic diffusion in disordered sys¬ 
tems, such as the many examples of subdiffusion recently 
observed in living cells [12-16, 38]. Notably, our formu¬ 
lation includes the possibility to model the simultaneous 
occurrence of subdiffusion (as well as any other types of 
motion) at the single particle level (Fig. 1) and EB (Fig. 
2 ), a feature observed in many experimental reports [12- 
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15]. This is in contrast with other nonergodic models, 
such as the CTRW, predicting a linear scaling of the 
time-averaged mean-squared displacement. Therefore, 
the data could not be satisfactorily interpreted by the 
CTRW alone and needed to include an additional source 
of subdiffusion together with CTRW models [14, 15]. 

In particular, we showed that our framework offers an 
interpretation of the data of Golding and Cox [12] on the 
basis of a fBm in a heterogeneous medium. The stochas¬ 
tic process (7) allows to capture both the subdiffusiv- 
ity in the time-averaged mean-squared displacement, the 
monotonic temporal growing of the p-variations test (as 
for the fBm), as well as the EB parameter value of the 
CTRW. Therefore, our model allows to reproduce all the 
features observed experimentally and thus solve the dis¬ 
agreement about the underlying stochastic process. 

Furthermore, we show that by introducing a nonsta¬ 
tionary random medium (19), our model can be extended 
to include the occurrence of aging, a feature often as¬ 
sociated to EB in living systems [14, 16]. As such, we 


consider that our general approach could contribute to 
investigate the occurrence of EB and anomalous diffu¬ 
sion in life sciences as well as many other fields, and help 
to elucidate their effects and implications. 
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FIG. 2. Plot of Eb{T) for the stochastic process (7) with 
P = 0.5 and H = 0.3 at various timelags A as a function of 
the measurement time T. Larger A produces an increase of 
Eb(T) at short time T. The Eb(T) values at large time T 
show are in agreement with the theoretical expectation (16) 
(dashed line). (Inset) Results of the the p-variation test with 
p — 2 for the stochastic process (7) with /3 = 0.5 and H = 0.3, 
showing the same trend as the pure fBm. 


FIG. 1. Time-averaged mean-square displacement 5^ as a 
function of the timelag A calculated for several trajectories 
(thin red lines) performing the fBm in a random medium, 
according to the ggBm (7) with P = H = 0.3 and T = 10^. 
Dashed line corresponds to the time and ensemble averaged 
mean-square displacement. Continuous thick line is a guide 
to the eye. 
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FIG. 3. Plot of the time and ensemble averaged mean-square 
displacement {S'^{T)} at various timelags A and as a fucntion 
of the measurement time T for the process (19) with (3 = H = 
0.3 and a = —0.3. The curves asymptotically show a power 
law decay T°‘ (dashed lines) demonstrating the presence of 
aging in the process. 
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FIG. 4. Contour plot of the aging exponent a as a func¬ 
tion of the exponents controlling the power law growth of the 
time- {2H) and ensemble-averaged {a + 2H) mean-squared 
displacement for the process (19). The continuous black line 
corresponds to the absence of aging {a = 0). Dashed green 
lines separate sub- and super-diffusive regions, characterized 
by exponents values < 1 and > 1, respectively. 



